The higher-dimensional Ablowitz-Ladik model: from (non-)integrability and solitary 
waves to surprising collapse properties and more exotic solutions 
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We propose a consideration of the properties of the two-dimensional Ablowitz-Ladik discretization 
of the ubiquitous nonlinear Schrodinger (NLS) model. We use singularity confinement techniques 
to suggest that the relevant discretization should not be integrable. More importantly, we identify 
the prototypical solitary waves of the model and examine their stability, illustrating the remarkable 
feature that near the continuum limit, this discretization leads to the absence of collapse and complete 
spectral wave stability, in stark contrast to the standard discretization of the NLS. We also briefly 
touch upon the three-dimensional case and generalizations of our considerations therein, and also 
present some more exotic solutions of the model, such as exact line solitons and discrete vortices. 
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Introduction. The nonlinear Schrodinger (NLS) equa- 
tion [J 0] is a prototypical dispersive nonlinear partial 
differential equation (PDE) that has been central for al- 
most four decades now to a variety of areas in mathemat- 
ical physics. The relevant fields of application vary from 
optics and propagation of the (envelope of) the electric 
field in optical fibers [3, IH , to the self-focusing and col- 
lapse of Langmuir waves in plasma physics [f| and from 
the behavior of deep water waves and freak waves in the 
ocean [f| @] to the mean-field dynamics of Bose-Einstein 
condensates in atomic physics [8, 9]. 

Much of the NLS literature and physical investigations 
have been centered around the one-dimensional (Id) set- 
ting, due not only to its mathematical and computational 
simplicity, but also in major part due to its complete inte- 
grability via the inverse scattering transform Yet, the 
higher dimensional investigations of the NLS equation 
have important elements of mathematics and physics in 
their own right, presenting the possibility for self- focusing 
and wave collapse [l| that has remarkable manifestations 
in some of the physical areas represented above. It is, 
thus, not surprising that recent fundamental investiga- 
tions have been focused on related issues, such as e.g. the 
observation of the self-similarly collapsing solitary wave 
of the two-dimensional (2d) NLS equation (the so-called 
Townes soliton) and on how to avoid the relevant 
collapse phenomena by means of temporal [ll|, [12] or 
spatial [131 ] variations of the nonlinearity, or by posing 
the problem on a lattice [Hj ]. 

In the present work, we revisit the higher dimensional 
NLS equation and consider its properties on a spatial 
lattice. This is an interesting approach in its own right, 
at least in part because many of the problems associ- 
ated with the optics (namely, optical waveguides 15] ) 



or atomic physics (namely, Bose-Einstein condensates - 
BECs- in optical lattices [16[) are inherently associated 
with such discrete models. Another relevant motivation 
is that of computation, since even computing with the 



continuum model proper requires posing the problem on 
a computational grid (with a finite, but small spacing). 
However, in our investigations herein, we will not use the 
"canonical" discrete form of the NLS equation, the so- 
called DNLS model 17, [3. Instead, we will consider 



an unconventional discretization that has yielded consid- 
erable insight (although, for a quite different reason, as 
illustrated below), namely the Ablowitz-Ladik (AL-NLS) 
model which in its two-dimensional form reads: 

^n,m — £ (^n+l,m "k ^n— l,m ~t~ ^n,'m+l ""T" ^n,m— 1 4ti Ujm ) 

+ ^ M [Un+l,m + Un-l,m + U n ,m+1 + Wn,m-l) (1) 

In particular, the fundamental difference from the DNLS 
here is that in addition to the centered-difference approx- 
imation of the Laplacian (with e = 1/Ax 2 being the 
coupling strength), a nearest-neighbor average is used 
to discretize the cubic nonlinearity of the continuum 



limit o~\u\ u (instead of a local term o~\u r , 



in the 



DNLS). The cubic nonlinearity physically represents the 
Kerr effect in optics (i.e., the dependence of the refractive 
index of the optical material on the light intensity) or the 
mean-field approximation of interatomic interactions as 
a nonlinear self-action in BECs. 

Our main findings and presentation are as follows. In 
the next section, we use the technique of singularity con- 
finement to illustrate that, by analogy to the continuum 
model, the AL-NLS model is unlikely to be integrable in 
two spatial dimensions. We then embark on a systematic 
analysis of the model's solitary waves and their stability 
and substantiate a surprising result, namely that the AL- 
NLS discretization possesses spectrally stable solitons for 
arbitrarily small lattice spacings Ax. We then consider 
the generalization of this result in the 3d setting. Lastly, 
we present some more exotic solutions of the AL-NLS 
model, such as the analytically exact (but unstable) line 
soliton and the x-shaped discrete vortex. 

Singularity Confinement. To determine whether the 
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AL-NLS is completely integrable, we use an integrability 
detector designed for difference or differential-difference 
equations, namely the singularity confinement (SC) cri- 
terion 1!J 2(J. The SC deals with the spontaneous ap- 
pearance of a singularity at some point in the lattice. The 
criterion is satisfied if the singular behavior is confined 
in a finite region of the lattice. In order to apply the SC 
criterion, we re-interpret the AL-NLS as an iteration for 
the real and imaginary parts of u = v + iw: 



W n +l,m = 



, m 
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If, at some point (jiq, mo), the denominator above is zero 
for some time t = t Q , then the iterates v n+ i m , w n+ i, m 
will be singular. We thus set 4e — a\u no ^ mo | 2 = (t—t )a(t) 
(for an arbitrary function a nonzero at t — t ), and let 
v n}m and w n>m for n — uq (with m ^ mo) and n — 
iiq — 1 to be regular and arbitrary. Computation of the 
iterates indicates that the singular behavior at t = t is 
not confined but rather propagates indefinitely for larger 
values of n. This result implies that the AL-NLS is not 
completely integrable. 

Solitary Waves and Stability. We start by recalling 
that the continuum analog of the 2d model is unsta- 
ble (due to an instability which is weaker than expo- 
nential (2l|) towards self-similar collapse [lj. The soli- 
tonic standing wave continuum solution (in the form 
u(x,y,t) = exp(iAt)v(x,y), taking Ax — > and a = — 1 
in Eq. (TT])) is the so-called Townes soliton [22J with a 
squared L 2 -norm (i.e., optical power or number of atoms 
in BEC) equal to P c w 11.7. This is the separatrix be- 
tween the collapse regime, arising for powers P > P c , 
and dispersion occurring for P < P c . 

In the discrete case of the AL-NLS, we also seek stand- 
ing waves in the form u n m (t) = exp(iAt)v n , m . Fixing 
a = — 1 (i.e., focusing nonlinearity) , we can seek such 
solutions either as a function of A or (equivalently, up 
to a simple rescaling) as a function of Ax. We present 
both cases in Fig. [1] The continuum limit is obtained as 
A or, respectively, Air — + 0. The former case has a par- 
ticularity with respect to stability (which is part of the 
reason as to why it is presented herein). In particular, 
the so-called Vakhitov-Kolokolov criterion suggests that 
for fundamental solutions [l], [13, EH ; the soliton stability 
is solely determined by the sign of the quantity dP/dA. 
In the AL-NLS case, interestingly (a direct generalization 
of the corresponding Id conservation law [2|), the power 
is defined as: 

c, / a , 
-In 1 \v 

4 V 4e 1 

such that it gives the continuum power in the limit. 
When, in our units, this quantity is positive, the soli- 
tary wave should be stable, while its negativity should 
imply exponential instability. 
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It is perhaps worthwhile to commence our comparison 
of the discrete solitary wave results to the corresponding 
continuum model through the intensely studied DNLS 
model 17, EH. There it is known (see also the right 
panel of Fig. [Tj for the relevant properties as a function 
of Ax) that as the continuum limit is approached, at the 
critical threshold of e = A (for spacings or frequencies 
below that), the solitary waves become exponentially un- 
stable. This dynamical instability, manifested the right 
panels of Fig. [2] for two different values of the spacing 
Ax, is connected to the collapse of the continuum model 
and consists of a "quasi-collapse" phenomenon, whereby 
all the power of the solution is "collected" on a single 
site; the discrete model with its discrete power conserva- 
tion and its spatial grid scale disallows a true self-similar 
collapse leading, in principle, to a finite time singularity. 
It should be noted that once unstable, the fundamental 
soliton of the DNLS model remains unstable throughout 
the parametric continuation to the continuum limit. The 
spectral picture of the linearization around the solitary 
wave (to examine the fate of small perturbations) reveals 
that in addition to a pair of eigenvalues always at the 
origin [due to the U(l) phase invariance of the model], 
there are 6 other pairs near the origin, approaching it, 
as the 2d continuum limit draws near. Of these, two 
identical pairs are associated with the translations along 
the x- and y-directions (whose corresponding invariances 
are restored in the limit) and which always remain stable 
(i.e., with A 2 < 0). On the contrary, the eigenvalue pair 
with A 2 > connected with the slope condition and with 
the quasi-collapse (focusing instability in the language of 
[23|]) is also connected with the additional symmetry of 
the 2d limit, namely the so-called pseudo-conformal in- 
variance [l[ , which allows the self-similar reshaping of the 
solution (and, hence, the access to the collapse dynamics 
of the continuum limit). 

One of the remarkable findings of the present work 
is that AL-NLS model are the fundamentally different 
spectral properties of its solitary wave. In particular, as 
can be seen in the top panel of Fig. [TJ the change of 
sign of dP/dA (and, accordingly, the instability) arises 
for 0.34 < A < 0.87 for Ax — 1, or, respectively, the in- 
stability emerges for 0.4825 < Ax < 0.758 (for A = 1.5). 
However, as the continuum limit is approached cither 
through A — > 0, or through Ax — > 0, the relevant eigen- 
value associated (in the limit) with the pseudo-conformal 
invariance has A 2 < 0. Furthermore, the squared eigen- 
values double pair associated with translations also ap- 
proaches zero from the negative side. This implies that 
arbitrarily close to the limit, the AL-NLS discretization 
offers an alternative of NLS free of collapse instabilities. 
This feature is also evident in the case of the dynamical 
evolution of Fig. [21 whose top left panel illustrates that 
for small values of Ax, the solitary waveform is not sub- 
ject to collapse (as it may be e.g. within the instability 
interval arising for higher A X, (IS illustrated in the figure 
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FIG. 1: (Left) Dependence of the AL-NLS fundamental soli- 
ton power P (Eq. (J2Jl) and of the square eigenvalues of the 
linearization A 2 on the wave frequency A (Aa; = 1). Vertical 
red lines bound the instability region of dP/dA < 0. Green 
dashed line denotes the (lower) edge of the phonon band. 
(Middle) Same but for the variation as a function of the grid 
spacing As: (A = 1.5). (Right) Same as the middle panel, but 
now for the regular DNLS discretization. 



for Ax = 0.6). 

Further consideration of this feature suggests that it is 
a particular trait of critical settings, which are at the very 
special (yet, physically realizable and important) separa- 
trix between subcritical settings where the solitary waves 
are dynamically stable and supercritical ones, where the 
waves are exponentially unstable. In this critical case, 
the linear spectrum possesses an additional zero eigen- 
value pair (associated with the pseudo-conformal invari- 
ance), which permits the reshaping of the solution under 
the action of the group of rescalings, and hence paves the 
way for the emergence of self-similar collapse. Discrete- 
ness can then shift this pair along the imaginary axis or 
along the real axis (as it is well documented to poten- 
tiall y d o with respect to translational eigenvalues also 
[H, llll)- The AL-NLS discretization turns out to be 
a prototypical example whereby the eigenvalue pair for- 
merly associated with the pseudo-conformal invariancc 
is perturbed in a stable way (moves along the imaginary 
axis of the spectral plane), upon discretization and hence, 
this model allows infinitesimally small spacings to give 
rise to collapse-free dynamics. 

Higher- dimensional Generalizations and More Com- 
plex Waveforms. The above discussion paves the way 
for understanding the approach to the continuum limit 
of the 3d solitary wave in the 3d AL-NLS generalization. 
The relevant power-frequency diagram is shown in Fig. 
[31 clearly illustrating through its negative slope the com- 
pletely unstable approach to the limit (since dP/dA < 0) 
In this case, the continuum limit has P « 18.82 [24j 
and the soliton is strongly (exponentially) unstable due 
to the supercritical collapse. Hence, the discretization 
has no choice but to respect the relevant limit and thus 
infinitesimally small grid spacings are unable to provide 




FIG. 2: Dynamical Evolution of the soliton solution max- 
imum for AL-NLS (left) and DNLS (right) for As; = 0.2 
(top) and 0.6 (bottom). The insets show the contour plot of 
|u n , m (100)| 2 . Notice the stability of the AL-NLS for Aa; = 0.2 
(all other cases lead to collapse) 
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FIG. 3: Power vs. frequency plot for the unstable dP/dA < 
3d AL-NLS case (Aa: = 1). 



stabilization irrespectively of the specific form of the dis- 
crete model. 

Lastly, we consider some more complex waveforms that 
can arise in the AL-NLS model. It is possible to obtain 
explicit solutions in the form of quasi-ld solitons e.g., 



^4sech(cm + f3m + xq) exp(iAt) 



(3) 



where the parameters satisfy A — ± y/ (A 2 + 8 Ae) / (4e) , 
and a = j3 = cosh - [(A + 4e)/(4e)]. Such a solution 
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FIG. 4: The top panels show the exact line soliton con- 
tour (left) and its dynamical instability evolution (right); 
A = 2s = 1. The bottom ones show the vortex square modu- 
lus and phase for A = 1.1 (left) and its power and linearization 
spectrum similar to Fig. [T] (right). 

is depicted in the top left panel of Fig. |H but the lin- 
earization around it illustrates that it is highly unsta- 
ble and its dynamics spontaneously lead to filamentation 
and the formation of localized solitary waves of the type 
considered previously, as seen in the top right panel of 
Fig. 2J On the other hand, there also exist more com- 
plex solutions, such as the discrete x-shaped vortices of 
the bottom left panel of Fig. [4] Such solutions exist in 
the DNLS equation, because of its anti-continuum limit 
e = [25j but do not persist in the continuum limit, 
hence their existence is not guaranteed in the AL-NLS 
model. Nevertheless, we find here that they exist and 
are quite robust, becoming unstable due to a real eigen- 
value pair for Ax < 1.25 (bottom right panel of Fig. Q|. 
Additional eigenvalue pairs emerge for Ax < 0.985. It is 
interesting to note that the corresponding x-shaped vor- 
tex in the DNLS case becomes unstable due to complex 
eigenvalue quartets for spacing values Ax < 1.178. 

Conclusions. In this study, we considered the 
Ablowitz-Ladik discretization of the NLS in higher di- 
mensional settings. We illustrated, via singularity con- 
finement, that the model is unlikely to be integrable, yet, 
due to the critical nature of the nonlinearity in 2d, it 
possesses some remarkable features, including robust and 
spectrally stable solitary waves even infinitesimally close 



to the continuum limit, explicit analytical solutions and 
more complex vortex waveforms. This study raises many 
fundamental issues. It would be relevant to obtain a sys- 
tematic understanding of how different NLS discretiza- 
tions in the critical case affect the pseudo-conformal in- 
variance (and whether there might conceivably exist one 
that preserves the symmetry). It would also be useful 
to examine how different types of discretizations affect 
other classes of critical models, such as the critical gen- 
eralized KdV . Lastly, it would be relevant to obtain 
a systematic classification (analogous to the one existing 
in DNLS of the solutions of the AL-NLS and their 
stability properties, in two- and higher dimensions. 
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